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Abstract 

Head-on collisions between two solitons in the pure CP 1 model on a flat torus are investigated 
via numerical simulations. The charge-two lumps, written out in terms of Weierstrass' elliptic 
p- function, are found to scatter at 90°. The phenomenon of singularity formation is also seen. 



1 Introduction 

The non-linear 0(3) sigma or CP 1 model in three-dimensional space-time is a rich industry 
of research, both for its condensed matter applications and as a simple field theory possessing 
topological solitons. The energy density associated with the solitons are lumps localised in space. 
The model also appears as a low dimensional analogue of non-abelian gauge field theories in (3+1) 
dimensions, an example being the increasingly popular Skyrme model of nuclear physics. In pure 
mathematics the CP 1 solitons are known as harmonic maps, by itself a long established area of 
research. 

The classical (2+0)-dimensional CP 1 model on the extended plane or Riemann sphere K2 U 
{00} ~ S2, where the solitons are harmonic maps S2 *— > S2 given by any complex holomorphic 
function, has been amply discussed in the literature jj], [|. The full time-dependent model [(2+1) 
dimensions] is not integrable, so numerical simulations are needed for studying its dynamics. Re- 
garding collisions, it is well known that the lumps on the extended plane scatter-off at 90° with 
respect to the initial direction of motion in the centre-of-mass frame M] . Due to the conformal in- 
variance of the planar model, the lumps are unstable in the sense that they can shrink indefinitely, 
leading to singularity formation in finite time Q. This instability is cured by supplementing the 
0(3) lagrangian with a Skyrme- like and a potential-like term 

More recently, attention has also been paid to the CP 1 model with square periodic boundary 
conditions, where the solitons are harmonic maps T2 1— > S2 (T2 a flat torus) given by elliptic 
functions. Physically, this approach looks more attractive than the one on 52 because the system 
is located in a finite volume from the outset. Mathematically, the model on T2 has the formal 
advantage of removing the problem confronted in the extended plane, whose non-compactness 
brings about difficulties to defining the metric on the moduli space of static soliton solutions. 

In reference Q| we tackled the toroidal model by expressing the soliton fields in terms of Weier- 
strass' er-function. We found several properties not seen in the usual planar theory, principally 
that there are no single-soliton solutions on the torus and, in the Skyrme version of the model, 
that there is no critical speed below which the skyrmions do not scatter at 90°. 
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Our anterior paper M resorted to the Weierstrass' elliptic p-funcion to numerically study the 
CP 1 model, both in its original and Skyrme versions, for systems with no initial speed. Novel 
features were unveiled in the topological charge two sector. Amongst them: • the appearance of 
four energy chunks instead of two for certain configurations initially situated in the diagonals of 
the lattice and • lump splitting in the Skyrme case for systems initially located in the central 
cross of the grid. 

Within the framework of the geodesic approximation, the elliptic p-function was employed as 
well by Speight to describe pure CP 1 lumps. His analysis predicted that the solitons may 
shrink and form singularities in finite time -as their siblings on S2- and that the solitons should 
also scatter at 90°. 

In the present paper we report the results for CP 1 lump scattering obtained by numerically 
simulating the initial-value problem given by the Weierstrass' p-function in the original, unmod- 
ified CP 1 model. Our results bear out the predictions referred to in the previous parragraph. 

The paper is arranged as follows: In the next section we lay out the CP 1 model with periodic 
boundary conditions. In section 3 the numerical procedure is explained and in section 4 the 
scattering results are discussed. We close the paper with some concluding remarks including 
suggestions for further research. 



2 The CP 1 model on the torus 

Our model is defined by the lagrangian density 



\d t W\ 2 -2\d z W\ 2 -2\d,W\ 2 

(l + \W\ 2 ) 2 1 ' 

where z = x + iy £ T?; z is the complex conjugate of z. The complex field W obeys the periodic 
boundary condition 

W[z+(m + in)L]=W(z), Vt, (2) 

where m,n = 0,1,2,... and L is the size of the torus. The static soliton solutions are elliptic 
functions which may be written as 

W = \p(z - a) +b, X,a,b € Z, (3) 

being p(z) the elliptic function of Weierstrass. Within a fundamental cell of vertices 

(0,0), (L,0), (L, L), (0,L), 

p is expandable as [g| 

p(z) = z - 2 + ^ z 2 + ^ z 4 + ...+^z 2j ' 2 + ..., 6 SR. (4) 

This function is of the second order (degree) and so (j|) represents solitons of topological charge 
2. As discussed in j^, (?J a distinctive feature of the toroidal model is the absence of analyt- 
ical single-soliton solutions. This can be understood by recalling that the simplest non-trivial 
elliptic function is of order 2. In the language of differential geometry, the degree of the har- 
monic maps M 1— > S2 must be greater than genus(M), M a compact and orientable Riemann 
surface. Since genus(T-2) = 1 there are no unit-degree maps on the torus. However, a periodic 
single-soliton ansatz was constructed in Q using the equation (^) below with k=1 and relaxing 
the accompanying selection rule. Note that when M = S2 we have the common model with 
standard boundary conditions which do possess solitons in all homotopy classes; this is because 
genus(S2) = 0. 

In reference || we computed the periodic solitons through 

3=1 3=1 3=1 
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with a subroutine that calculates <j(z) numerically. Although a is pseudo-elliptic, the above ratio 
subject to the summation of its zeros (a,j) being equal to the summation of its poles (bj) renders 
W elliptic. In this paper we utilise a similar subroutine and then reckon p(z) using the formula 



p(z) = -^ln[a(z)]. (6) 



The Laurent expansion for a reads 



«r(*) = °i e & (7) 

3=0 



Assisted by the useful properties of p(z) 



PO) = P(~z), p( iz ) = -p(z)j p( z ) = p( z ), (8) 

one readily deduces that p(z) is real on the boundary and central cross of the fundamental cell, 
and purely imaginary on the borders JhJ. 

The static (or potential) energy density E associated with the harmonic map W may be 
read-off from the lagrangian (|l|). Taking into account the identity 

[^LlAf = Mz - a)[p(z - a) 2 - p a (L/2)], (9) 

we have 

E = 8|A| [l + \Xp(z-a) + b\^ • (10) 

3 Basic numerical procedure 

We treat configurations of the form (^) as the initial conditions for our time evolution, studied nu- 
merically. The time derivative of W is calculated from the Lorentz-boosted field. Our simulations 
run in the ^-formulation of the model, whose field equation 

= [-(dttf + (a^) 2 + {a y ^f]~$ + dU + (ii) 

follows from the lagrangian density ([!]) with the help of 

W =l^-. (12) 

The real scalar field ~~$ — ((pi, 4>2, $3) satisfies ~~$ -~$ = 1. Inverting formula ([[2]) entails 

^_ ( w + w . -W + W \W\ 2 -1 

9 l |W| 2 + l' \W\ 2 + 1' \W\ 2 + 1>' ( ' 

We compute the series (Q) up to the fifth term, the coefficients Cj being in our case negligibly 
small for j > 6. We employ the fourth-order Runge-Kutta method and approximate the spatial 
derivatives by finite differences. The laplacian is evaluated using the standard nine-point formula 
and, to further check our results, a 13-point recipe is also used. Our results show unsensitiveness 
to either technique, thus confirming the reliability of our results. The discrete model evolves on a 
200 x 200 periodic lattice (n x = n y = 200) with spatial and time steps 8x—Sy— 0.02 and 5t=0.005, 
respectively. The size of our fundamental, toroidal network is L — n x x 5x — 4. 

Unavoidable round-off errors gradually shift the fields away from the constraint <j> . <j> =1. So 



we rescale <j> — > <f> /y 4> ■ 4> every few iterations. Each time, just before the rescaling operation, 
we evaluate the quantity /1 = 0.0—1 at each lattice point. Treating the maximum of the 
absolute value of /1 as a measure of the numerical errors, we find that max|/x w 10~ 10 . This 
magnitude is useful as a guide to determine how reliable a given numerical result is. Usage of 
an unsound numerical procedure in the Runge-Kutta evolution shows itself as a rapid growth of 
max|/i|; this also occurs, for instance, when the unstable energy lumps become too spiky. 
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4 Results 



The initial conditions are then given by equation (|3|) 

W = A p(z - a) + b, 

where A is related to the size of the solitons, 6 determines their mutual separation and a merely 
shifts the solution on the torus. Calling the initial speed v, we boost the above field and take its 
time derivative: 

W -> W(t) = Xp(z-a)+b(l-vt), 

d t W(t) = ~bv, { ' 

thus completely defining our initial-value problem. Without loss of generality we may take the 
values of these parameters according to numerical convenience. Let us set 

A = l, a = (2.015,2.015), 6=1 

throughout. As depicted in figure [l], the total energy density corresponding to our soliton field 
has the form of two lumps sitting on the central cross of the cell, symmetrically around its centre. 

Consider the situation where the lumps are sent towards one another with a relative initial 
speed of v = 0.3. We observe that the solitons gradually expand as they approach each other; 
they collide at the centre of the nett and coalesce into a ringish structure, where the solitons 
are no longer distinguishable. Here they attain maximum expansion, i.e., the peak of the total 
energy density (E max ) reaches a minimum value. After this process the lumps get narrower and 
narrower as they re-emerge at right angles to their initial line of approach. They keep shrinking 
while moving away from the centre, in opposite senses. Some time later they become so spiky 
that the numerical code breaks down. This is the well-documented instability of the 0(3) model 
in two dimensions, where the theory in conformally invariant. The foregoing events are illustrated 
in figure ^. 

The evolution of E max is shown in the upper half of figure ^ for two values of v. The curves 
are qualitatively alike, the life of the system being longer the smaller its initial velocity is. 

The kinetic energy remains very small all along except when the instability takes over. This 
can be appreciated from figure ^| (bottom), where the maximum value of the kinetic energy density 
(Km ax ) is plotted versus time. 

Solitons with no initial velocity remain motionless with the passing of time (see below, however, 
the particular case when 6 = 0). They too exhibit shrinking and singularity formation, if at a 
slower rate than when v 7^ 0; the evolution of the total energy density qualitatively resembles that 
of figure ^ (top). As for the kinetic energy density, it stays very small, at around w 10~ 7 during 
the whole simulation (except when the breadth of the lumps is comparable to the lattice spacing 
which leads to the collapse of the numerics). 

Particularly noteworthy is the case 6 = 0. The soliton field and its energy distribution are, 
respectively from equations (||) and (|Io|), 

W = \p(z — a), 

E _ g| A |2 |p( Z -a)||p 2 ( Z -a)-^ 2 (L/2)| (15) 

The global maxima of such energy density are located along the diagonals of the basic cell || , 
i.e., where p[z) is purely imaginary. Given that E in (|l^) is invariant under p —> — p, it follows 
from the evenness of p that E has at least for peaks on the diagonals of the fundamental cell. In 
fact, it possesses four peaks (see figure ^). This situation in the topological index- two class is a 
remarkable feature of p-solitons, going beyond the two-lump and annular structures found in the 
CP 1 on S 2 . 

Now, the scattering problem when 6 = is especial: In order to zoom the lumps towards each 
other one must boost W according to equation (p^), i.e., set 6 — » 6(1 — vt). But this is only 
possible if 6 0. And boosting the parameter a in W is of no avail, for it simply shifts the system 
as a whole. Let us then consider systems of the form (|l^) started off from rest. It turns out that 
such configurations shrink quite slowly with time, as illustrated in the upper graph of figure ||. 
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From the bottom-right plot of figure |5| we observe that the solitons move away from the centre 
and back in along the diagonals of the grid in breather-like fashion. The kinetic waves emitted in 
the process evolve according to the bottom-left side of illustration number ^. 

Given similar set-ups, CP 1 lumps with periodic boundary conditions will usually shrink at 
a slower rate than their relatives with standard boundary conditions. This is because solitons 
in compact domains as T2 always coalesce somehow, and these are configurations with E max 
smaller than when the solitons do not overlap -a situation achivable to a good approximation in 
non-compact spaces like the extended plane. States with b = 0, which correspond to coincident 
structures that spread themselves all over the cell (see picture^), have a small E max which evolves 
quasi-stably as depicted in figure ^. Note that other than through the parameter b, we can always 
bring the solitons to greater overlapping either by increasing the lump size (A) or by reducing the 
torus size (L). This is be useful if we want to run a collision experiment with the same value of b 
but different initial values of E max . 

We point out that in reference Q our configurations for 6 = exhibited scattering in the 
Skyrme version of the model, even though the initial speed was zero. We define our Skyrme 
model on the torus by [compare with the lagrangian (ill)] 



\dtW\ 2 - 2\d z W\ 2 - 2\dzWf_ 



(1 + \W\ 2 ) 2 



|cW 
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(16) 

In this scheme the four energy peaks move and collide along the diagonals and scatter at right 
angles. Let us emphasise that the extra term in the lagrangian ([Hj) stabilises the CP 1 lumps, 
i.e., prevents the formation of singularities. 

As pointed out in the introduction, our main results for p-solitons (singularity formation and 
scattering at 90°) conform with prognostications made in the geodesic approximation treatment of 
the model There it was proved that the moduli space of static 2- lump solutions is geodesically 
incomplete and has finite diameter, leading to infinite lump shrinkage. It was also proved in 
[^] that some mathematical constraints oblige the lumps to lie either on the central cross, the 
boundary or the diagonals of the lattice, implying that only scattering at right angles should 



5 Conclusions 

With the help of numerical simulations we have investigated some scattering properties of the 
(2+1) dimensional nonlinear 0(3) or CP 1 model with periodic boundary conditions, expressing 
the soliton solutions through Weierstrass' elliptic p-function. Limiting ourselves to the pure 
version of the model, we have observed that the lumps scatter off forming 90° with the initial 
direction of motion in the centre-of-mass frame. During this process the lumps grow spiker as 
time elapses, eventually breaking down the numerical code. Lack of stability shows as well for 
lumps with zero initial speed, but at a slower rate. The 0(3) instability may be understood from 
the symmetry of the model under dilation transformations. The above properties are basically the 
same as those known from the familiar model on the Riemann sphere. But some other phenomena 
turn out to be quite disimilar on T2, e.g., the absence of unit-charge solitons and the presence of 
four lumps rather than two in the charge-two topological class. The non-existence of single-soliton 
fields on T2 is dictated by a general property stating that there are no elliptic functions of order 
one. The four peaks in the charge-two sector stem from a notable symmetry property specific to 
the p-energy density when b = [see equation ([l5])] . A configuration of this type is not observed 
when the fields are constructed out of the er-function. 

A natural extension of the present work is the study of p-lumps with non-zero initial velocity 
in the Skyrme format. Skyrmions with such initial conditions have already shown a thought- 
provoking scattering and splitting pattern under numerical simulations. Consequently, research 
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Figure 1: Three dimensional and contour pictures corresponding to the total energy density of the 
solitons at t = 0. The initial speed is v — 0.3. 
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Emax=1 6.4; t=3 





Emax=64.5;t=10 



Figure 2: Continued. The lumps scatter at 90° to the original direction of motion. 
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Lum ps at the initial tim e forthe case b=0 



Emax=10.94 




Figure 4: Lumps at t = corresponding to the particular case W = p{z — a). The presence of four 
lumps rather than two in the charge two sector is an alluring characteristic of p-solitons. 
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Graphs for the b=0 case 
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Figure 5: Plots corresponding to W = p(z — a). We observe a quasi-stable configuration where the 
lumps oscillate along the diagonals of the network in breatherdike manner. The initial velocity is zero. 
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on their evolution under boosting seems compelling. Further investigations may involve defining 
the solitons in terms of other elliptic functions as those of Jacobi, and the pseudo-elliptic 9- 
functions. A classification of solitonic properties on the torus can only be made after thorough 
consideration of the above-mentioned functions. Let us remind that on S2 no new traits arise 
from casting the soliton configurations in different ways. In the topological charge two class, for 
example, the following fields behave qualitatively alike on the Riemann sphere: 

(2 — cj(z — d) 

It is a non-trivial problem to understand the mechanisms underlying the complicated CP 1 dy- 
namics. In the quest for such understanding, a study of solitonic behaviour on compact manifolds, 
especially T 2 , is certainly quite an appealing program to pursue. 



12 



Acknowledgement 

The author is very grateful to WJ Zakrzewski for helpful conversations and to the 
Mathematical Sciences Dept. of the University of Durham, England, where part of the 
present work was carried out. The financial support of La Universidad del Zulia through 
CONDES 218-99 is acknowledged with big thanks. 

References 

[1] Eichcnhcrr H. (1976) Nucl. Phys. B146, 215 

[2] Pcrelomov A. M. (1981) Physica 4D, (1981) 1 

[3] Zakrzewski W. J. (1991) Nonlinearity 4, 429 

[4] Leese R. A., Peyrard M. and Zakrzewski W. J. (1990) Nonlinearity 3, 387 

[5] Cova R. J. (1995) Hclv. Phys. Acta 8, 282 

[6] Cova R. J. and Zakrzewski W. J. (1997) Nonlinearity 10, 1305 

[7] Cova R. J. (1998) Hclv. Phys. Acta 71, 675 

[8] Speight J. M. (1998) Comm. Math. Phys 194, 513 

[9] Goursat E. (1916) Functions of a complex variable, Dover Publications 

[10] Lawden D. F. (1989) Elliptic functions and applications, Springer Verlag USA 



13 



